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^ ! Abstract 

o 

f— ^ , We investigate the consequences of Fierz transformations acting upon the contact in- 

teractions for nucleon fields occurring in relativistic point coupling models in Hartree 
approximation, which yield the same models but in Hartree-Fock approximation in- 
stead. Identical nuclear ground state observables are calculated in the two approx- 
^ I imations, but the magnitudes of the coupling constants are different. We find for 

G ' model studies of four-fermion interactions occurring in two existing relativistic point 

coupling phenomenologies that whereas in Hartree the isovector-scalar strength 0^5, 
k><( \ corresponding to (5-meson exchange, is unnaturally small, indicating a possible new 

^H ' symmetry, in Hartree-Fock it is instead comparable to the isovector- vector strength 



axv corresponding to p-meson exchange, but the sum of the two isovector coupling 
constants appears to be preserved in both approaches. Furthermore, in Hartree-Fock 
approximation, both QCD-scaled isovector coupling constants are natural (dimen- 
sionless and of order 1) whereas in Hartree approximation only that of the isovector- 
vector channel is natural. This indicates that it is not necessary to search for a new 
symmetry and, moreover, that the role of the (5~meson should be reexamined. This 
work presents the first comparisons of naturalized coupling constants coming from 
relativistic Hartree and relativistic Hartree-Fock solutions to the same Lagrangian. 
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1 Introduction 



Relativistic mean field (RMF) models are quite successful in describing ground 
state properties of finite nuclei and nuclear matter properties. Such models 
describe the nucleus as a system of Dirac nucleons that interact in a relativistic 
covariant manner via mean meson fields [1,2,3,4] or via mean nucleon fields 
[5,6,7,8]. The meson fields are of finite range (FR) due to the meson exchange 
whereas the nucleon fields are of zero range (contact interactions or point 
couplings PC) together with derivative terms that simulate the finite range 
meson exchanges. A common element to the calculations referenced above is 
that they have all been performed in relativistic Hartree approximation. 

The RMF-FR studies to date have generally considered three explicit meson 
fields. These are the isoscalar-scalar field due to exchange of the a meson, the 
isoscalar-vector field due to u meson exchange, and the isovector-vector field 
due to p meson exchange. The isovector-scalar field due to 6 meson exchange 
has generally not been included because its contribution to the nuclear force 
from one-boson exchange is considered weak [10], given a relatively large mass 
of 983 MeV and a relatively small (but not well determined) coupling con- 
stant. It is, however, included in the RMF-PC studies of Refs. [5,6,8], where 
it is found that its contribution is very small. Furthermore, in a study [11] of 
the naturalness of the set of coupling coupling constants from Ref. [5] it was 
discovered that the isovector-scalar coupling constant is unnaturally small. 
This would presuppose a symmetry to preserve its small value. 

Thus, in relativistic Hartree approximation the isovector-scalar channel may 
be neglected, and (perhaps) a symmetry may be identified to preserve the 
small value of its coupling constant. However, we have not found any such 
symmetry. Therefore, we instead examine our calculational approach, the rel- 
ativistic Hartree approximation, and ask what happens to the magnitudes of 
the coupling constants in relativistic Hartree-Fock approximation where both 
the direct and exchange terms explicitly appear? 

We believe that the investigation of exchange terms in an effective field the- 
ory for nucleons is meaningful even if one considers the RMF model as an 
approximation to the exact density functional in the spirit of the Hohenberg- 
Kohn theorem and Kohn-Sham theory [9], where exchange effects should be 
absorbed in the various coupling constants. This is because the Hartree-Fock 
theory can be viewed as a Kohn-Sham formalism with exact treatment of ex- 
change (see Ref. [13], for example), which is then a different representation 
from the Hartree representation. It also has the correct one-particle limit and 
is a self-interaction free theory [12,13]. This is important for odd systems where 
the odd particle feels its own potential if exchange is ignored. The original idea 
to perform relativistic Hartree-Fock calculations by using contact interactions 



is due to Ref. [14] . 

In Sec. II we present the simplest possible Lagrangian containing scalar and 
vector fields of both isoscalar and isovector character, and we relate the cou- 
pling constants for this Lagrangian in relativistic Hartree approximation to the 
corresponding coupling constants in relativistic Hartree-Fock approximation. 
We apply our results in Sec. Ill to study the four-fermion contact interactions 
occurring in two existing realistic point coupling models determined (phe- 
nomenologically) in relativistic Hartree approximation. We then address the 
question of naturalness of the Hartree and Hartree-Fock coupling constants 
from these two models in Sec. IV. Our conclusions are given in Sec. V. 



2 Four-Fermion Relativistic Point Coupling with Exchange 



For comparing Hartree and Hartree-Fock representations we consider two- 
body contact interactions (four-fermion point couplings) in the mean field 
and no sea approximations, applied to the ground states of even-even nuclei: 



1 - 1 - 

- -^arsii^rij) ■ {^fip) - -aryl^T^^V') ■ {^l^fip) (1) 

where ip is the nucleon field and r is the isospin matrix. We wish to compare 
the same model ansatz in two different many-body approximations. Thus, our 
model space does not explicitly include pions because the pion field vanishes 
in the Hartree approximation, but contributes via its exchange terms in the 
Hartree-Fock approximation, which would then yield two different models. 
Implicitly the effects of the pion are nevertheless included because our coupling 
constants are determined by measured observables. Accordingly, we regard this 
work as a model study and do not construct a complete Hartree-Fock model. 

Taking the normal ordered expectation value of £ in a Slater determinant 1$) 
leads to the well-known direct and exchange terms that are to be solved in 
Hartree-Fock approximation, namely. 



($1 : Cnv ■ \^) = --asPs^ - -^aypv'^ - -^arsPrs^ - -aTvPTv"^ 
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+ 2^sPSeJ + 2^vpVeJ + -^arspTSeJ + ^(^TVpTVeJ (2) 



Here, we refer to the {a} as coupling constants, ps and pv denote the isoscalar 



scalar density and the time component of the isoscalar vector density, respec- 
tively, and pts and ptv denote the corresponding isovector densities. The 
squares of the isoscalar scalar and vector exchange densities, ps^^ and pv^^, 
and those of the isovector scalar and vector exchange densities, Pts^^ ^^^ 
PTVea:^ are given by [a and b are nucleon states]: 



PsJ = T.i^ai^b){Ma) (3) 

a,b 

PvJ = T.i^alt.A){^bY^a) (4) 

a,b 

Pts, J = ^(V^aT^b) ■ {■ipbTiIJa) (5) 

a,b 

Ptv J = ^^(V^aTM^V^fe) " {Al^rifja) (6) 

a,b 



We then apply Fierz transformations [15,16] in Dirac-iso space to the four 
exchange densities. This transformation expresses a product of nondiagonal 
matrix elements of Dirac F-matrices as an expansion into products of diagonal 
matrix elements, such as 
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{ijaTitjJb){ijbTjtpa) = J2 Ckl{^a^kiJa){^b^l^b) (7) 

k,l=l 



where Fj stands for one of the sixteen Dirac matrices {1,7^,75,757^,0-^1,} 
constituting a linearly independent basis in the space of complex 4x4 matrices, 
which may be coupled or uncoupled to isospin matrices r. Applying Eq. (7) 
to Eqs. (3-6), one immediately sees that all terms containing 75 and 757^ 
vanish because ipa and ipb are nucleon fields with good parity. We now make 
the (reasonable) approximation that the tensor (cr^j,) and iso-tensor (rcr^j,) 
contributions are quite small and can be neglected [7] . Reordering the resulting 
terms leads to a Lagrangian that is formally identical to Eq. (2) without the 
four exchange terms, but with newly defined coupling constants instead: 
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C~ = --asps^ - -OivPv'^ - -OiTSPTs'^ - -OiTVPTv' 



where the newly defined coupling constants are given by 



Ois = -^Ois - -ay - -axs - -jOirv 

15 3 3 

av = -o«s + Ttty - -^OLTs + jCiTv 

_ 1 1 9 1 

c^TS = ~'^'^s - -^av + -ars + -^o^tv 

1113 



This result already shows that in the Hartree-Fock approach, due to the ex- 
change effect, all original terms contribute to all channels of the effective inter- 
action. And the formal structure of Eq. (8) is identical to the Hartree approx- 
imation for the same model with redefined coupling constants. However, this 
Lagrangian, when considering all terms arising from the Fierz transformations, 
is a self-interaction free theory. The inverse solution of Eq. (9) is 
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Given the above results, if one determines the coupling constants of a rela- 
tivistic point coupling Lagrangian in Hartree approximation, then the set of 
coupling constants {a} in Eqs. (8) and (9) has been determined. Use of these 
coupling constants in Eq. (10) then yields the original Hartree-Fock coupling 
constants {a}. The two sets of coupling constants yield identical predictions of 
the nuclear ground state observables, but their magnitudes and physical inter- 
pretation are different because the former set implicitly accounts for exchange 
processes whereas the latter set explicitly accounts for exchange processes. 



Exchange Effects in Relativistic Point Coupling Models Deter- 
mined from Measured Observables in Hartree Approximation 



We examine coupling constants occurring in two realistic relativistic point 
coupling models that have been determined in Hartree approximation. The 
four terms of Eq. (1) are included, but higher order terms (six- and eight- 
fermion point couplings) and derivative terms are included as well. Here we 
focus our attention on the four four-fermion point couplings alone. 



Table 1 

Four-Fermion Relativistic Hartree {a} and Hartree-Fock {a} Coupling Constants 

in Two Realistic Lagrangians (PC-LA and PC-F4) [10"^ MeV"^]. 



Force 


as 


av 


ars 


drv 


as 


av 


ars 


arv 


PC-LA 
PC-F4 


-4.508 
-3.834 


3.427 
2.594 


7.421x10-3 
-5.924 xlO-2 


3.257 xlO-i 
3.937x10-1 


-5.712 
-4.608 


4.244 

2.872 


2.286 
1.345 


-2.314 
-1.425 



The two models are: PC-LA containing 9 coupling constants that appears 
in Ref. [5] published in 1992, and PC-F4 containing 11 coupling constants 
that appears in Ref. [8] published in 2002. As explained above, the four four- 
fermion coupling constants appearing in these tables are taken as the set {a} 
in Eqs. (8) and (9). The inverse solution Eq. (10) then yields the Hartree-Fock 
coupling constants {a}. We show both the Hartree {a} and Hartree-Fock {a} 
coupling constants in Table 1. 

Comparing the relativistic Hartree coupling constants of the models PC-LA 
and PC-Fl with the relativistic Hartree-Fock coupling constants one observes 
the following: (a) whereas the exponents in the Hartree coupling constants 
range from -7 to -4, those of the Hartree-Fock are all -4; (b) the Hartree-Fock 
isovector-scalar coupling constant is much larger than its Hartree counterpart, 
and has changed sign in PC-F4; and (c) the Hartree-Fock isovector-vector cou- 
pling constant has changed sign and its absolute magnitude has also increased 
in comparison to its Hartree counterpart. 

The Hartree-Fock isovector coupling constants have a larger role than those 
of Hartree. In fact, the four Hartree-Fock coupling constants are of the same 
order of magnitude and, furthermore, the magnitudes of the two isovector 
coupling constants are roughly equal. But none of this is true in the Hartree 
case. The Hartree-Fock coupling constants from both models have the same 
signs (not true for the Hartree case). In addition, it appears from the two 
models that the sum of the isovector coupling constants is better determined 
by the ground state observables than are the individual values, as was also 
learned in Ref. [81. 



Finally, we can ask what are the errors (uncorrelated and correlated) in the 
determination of the sets of coupling constants in the two calculations? We 
find that (a) the uncorrelated and correlated errors in the isovector coupling 
constants are significantly diminished in Hartree-Fock approximation (some 
by roughly two orders of magnitude), and (b) in this approximation, all four 
of the coupling constants are equally well determined, unlike in the Hartree 
approximation. This result is a consequence of the fact that in Hartree-Fock 
approximation the four four-fermion coupling constants contribute in each of 
the four channels, due to the explicit treatment of exchange processes. 



Therefore, the magnitudes of the four four-fermion Hartree-Fock couphng con- 
stants, in the two point-couphng models studied, are well determined and they 
are comparable. This immediately brings to mind the question of the natural- 
ness of these coupling constants to which we now turn our attention. 



4 The Quest for Naturalness 



The naturalness of the coupling constants relates to the question as to whether 
QCD scaling and chiral symmetry apply to finite nuclei. In 1990, Weinberg [17] 
showed that Lagrangians with (broken) chiral symmetry predict the suppres- 
sion of N-body forces. He accomplished this by constructing the most general 
possible chiral Lagrangian involving pions and low-energy nucleons as an infi- 
nite series of allowed derivative and contact interaction terms and then using 
QCD energy (mass) scales and dimensional power counting to categorize the 
terms of the series. This led to a systematic suppression of the N-body forces. 

We use the scaling procedure of Manohar and Georgi [18] but without pion 
fields. Explicit pionic degrees of freedom are absent in RMF Hartree theory, 
but can be present in RMF Hartree-Fock theory where Eq. (11) then also 
contains the pion field and pion mass as in Eq. (1) of Ref. [11]. The scaled 
generic Lagrangian term of the (physical) series is, without pions, 

C Cin ^^ — /^A (11) 



where -0 is a nucleon field, fj^ is the pion decay constant, 92.5 MeV, A = 770 
MeV is the QCD large- mass scale taken as the p meson mass, and {d'^) signifies 
a derivative. Dirac matrices and isospin operators (we use t here rather than 
t) have been ignored. Chiral symmetry demands [19] 



A = / + n-2>0 (12) 

such that the series contains only positive powers of (1/A). If the theory is 
natural [18], the Lagrangian should lead to dimensionless coefficients c^ of 
order unity. Our more stringent definition [8] is that a set of QCD-scaled 
coupling constants is natural if their absolute values are distributed about the 
value 1 and the ratio of the absolute maximum value to the absolute minimum 
value is less than 10. Thus, all information on scales ultimately resides in the 
c/„. If they are natural, QCD scaling works. 

Applying Eq. (11) to the dimensioned relativistic Hartree and Hartree-Fock 
couphng constants of Table 1, we obtain the corresponding sets of QCD- 



Table 2 

Relativistic Hartree {a} and Relativistic Hartree-Fock {a} Naturalized Coupling 
Constants {q„} for the Four-Fermion Point Couplings in Two Realistic Lagrangians 
(PC-LA and PC-F4). 



Coup. Const. 


Q„(PC - LA) 


Q„(PC-F4) 


as 


-L928 


-1.641 


OLy 


L466 


1.109 


ars 


0.013 


-0.101 


OITV 


0.557 


0.674 


# natural 


3 


3 


rnax / min 


152. 


16.2 


as 


-2.443 


-1.971 


ay 


1.815 


1.229 


ars 


3.912 


2.301 


arv 


-3.958 


-2.438 


# natural 


4 


4 


max /min 


2.18 


1.98 



scaled coupling constants listed in Table 2. This table shows that, whereas the 
Hartree coupling constants are not natural, the Hartree-Fock coupling con- 
stants are natural. The culprit is the isovector-scalar channel (corresponding 
to 6 meson exchange) in the Hartree approximation. Additional studies with 
toy Lagrangians consisting only of four-fermion interactions, not presented 
herein, show that naturalness is recovered in the Hartree-Fock representation 
from Hartree solutions where either of the isovector coupling constants is un- 
usually small. The culprit can apparently be either of the isovector coupling 
constants provided their sum remains approximately constant. The exchange 
process, however, dominates in both isovector channels to the extent that 
naturalness is recovered in the Hartree-Fock approximation in either case. 



The explicit inclusion of pion terms in the Hartree-Fock Lagrangian should 
not affect these results. This is because the generic chiral Lagrangian of Eqs. 
(11) and (12) has an important property: refining the model by adding new 
terms, such as pions, will change all of the naturalized coupling constants, 
but naturalness will still apply, that is, naturalness is largely independent of 
the details, such as adding pions, provided the physics is introduced via the 
measured observables in the framework of these equations. 



5 Conclusions 



We have extracted relativistic Hartree-Fock coupling constants from the cou- 
pling constants for relativistic Hartree calculations by use of Fierz relations 
and contact interactions. Identical observables are calculated with these two 
approximations, but the coupling constants and their physical interpretation 
are different because, whereas in the Hartree approximation they implicitly 
account for the exchange processes, in the Hartree-Fock approximation they 
explicitly account for the exchange processes. 

We have learned three things. First, viewing Hartree and Hartree-Fock as two 
different representations we have shown that the smallness of the isovector- 
scalar coupling constant axs Kiay be an artifact of the Hartree representation 
rather than the signature of a new symmetry and, therefore, the search for the 
symmetry is not needed. Second, while several relativistic Hartree studies have 
concluded that the observables determine a single isovector coupling constant 
(that of the p-meson), which is natural, we have learned that it is the sum 
of two isovector coupling constants that is preserved in relativistic Hartree 
and that individually these are both natural and their sum is also preserved 
in relativistic Hartree-Fock. Given that the (5-meson is now playing a role in 
recent asymmetric nuclear matter studies [20] it is relevant that its coupling 
constant is natural. Finally, the Hartree-Fock approximation may constitute a 
physically more realistic framework for power counting and QCD scaling than 
the Hartree approximation. 

Our conclusions are subject to the assumptions that: (a) we are able to neglect 
tensor contributions, (b) our results for four-fermion contact interactions will 
not be strongly affected by Fierz transformations on the remaining higher- 
order and derivative terms, and (c) explicit inclusion of pion interactions will 
change all of the coupling constants while retaining their naturalness. 
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